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Antiferromagnetism of the Hubbard Model 
on a Layered Honeycomb Lattice 
— Is MgB2 a Nearly-Antiferromagnetic Metal? — 
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As a model for the B 2pz bands of MgB2 , a half-filled Hubbard model on a honeycomb lattice 
is studied, and a possibility of an antiferromagnetism (AF) at the ground state is examined 
by Monte Carlo calculations as well as a random phase approximation method. On a single- 
layer honeycomb lattice, a paramagnetic to AF phase transition occurs at a finite value of the 
repulsion energy Uc/t ~ 3.6, due to the loss of the density of states at the fermi level. On a 
layered honeycomb lattice, however, inter-layer hopping term creates a perfect nesting of the 
fermi surface, and leads the AF with Uc = 0. From these facts, we discuss that MgB2 is possibly 
a nearly- AF metal, where the AF of the 2pz bands is destroyed by carrier doping. 
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§1. Introduction 

The recent discovery of MgB2j3^ an intermetaUic su- 
perconductor with high superconducting transition tem- 
perature (Tc), initiated a large number of activities both 
in experimental and theoretical researches-. Several ex- 
periments, including NMIl,® specific heatjj^ as well as 
photoemission spectrumo' measurements suggest that 
this material showSjJiGS^like s-wave superconductivity. 
Band calculationsEloQ'Q'B' indicate a strong interaction 
between electrons and high-frequency phonon modes ex- 
ists, which has been discussed as a possible origin of high 
Tc superconductivity in this material. It has been also 
argued that roles of strong Coulomb intcjjactiDns may 
directly or indirectly give rise to high rc.E3'llil'E3 

The crystal of MgB2 consists of alternative stackings 
of honeycomb-lattice B2 layers and triangular-lattice Mg 
layers, as is in AIB2. The band calculations show that 
Mg are almost fully ionized, and the bands at the fermi 
level have mostly B 2p characters. There exist four fermi 
surfaces, two of them in the 2pz bands and the rest in 
the 2px,y bands. The 2px.y bands are derived from the 
intra-layer a bonding orbitals, and have a quasi two- 
dimensional hole character. The 2pz bands are derived 
from the intra-layer tt bonding orbitals which also have 
inter-layer couplings between adjacent atomic orbitals in 
the c-axis direction. 

The -tight-binding fits to the band structure calcula- 
tionaj'O show that the 2pz bands are well reproduced 
by a nearest neighbor (N.N.) hopping model on a layered 



honeycomb lattice of B sites in the formulae 

e±(k) = ±2tj;yjxy{k) " 2t ^ cos(fc^), (1) 
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IxyCk) = y 4 + 2 + cos(-fcx) cos( — fcj,), (2) 

where txy {tz) denotes intra-layer (inter-layer) N.N. elec- 
tron hopping energy, while jxy is the structure factor 
of the intra-layer electron hopping to the N.N. sites on 
a honeycomb lattice. A lattice structure of the tight- 
binding model is illustrated in Fig. 1. 




Fig. 1. Boron sites in the layered honeycomb lattice structure 
of MgB2. Tight-binding fits show that a N.N. hopping model 
between boron sites reproduces the 2pz bands. Characters A 
and B denote sublattice indices. 



In this Letter, we investigate an aspect of the 2pz 
bands in MgB2 which comes from its characteristic band 
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structure. We consider a Hubbard model on a layered 
honeycomb lattice. Using a random phase approxima- 
tion (RPA) as well as quantum Monte Carlo (MC) meth- 
ods, antiferromagnetism (AF) at the ground state is ex- 
amined. We discuss that inter-layer electron hoppings 
create a perfect fermi-surface nesting in the 2^^ bands. 
With infinitesimal values of the Coulomb repulsion, the 
system is unstable against an AF order at half filling. 

First, let us consider a half-filled Hubbard model on a 
general bipartite lattice with N.N. hoppings. This model 
has a particle-hole symmetry, and as a result, a non- 
interacting electron dispersion has a relation 



e{k + Q) = -e{k), 



(3) 



where Q is the AF wavevector defined by the bipartite 
structure of the lattice, e.g. Q = (7r,7r) for a square 
lattice. To be more precise, the relation holds when 
a electron on a A-sublattice site hops only to sites on 
B-sublattice, and vice versa. A particle-hole symmetry 
gives the chemical potential /i = at half filling. On the 
fermi surface, the relation e(fc) — e{k + Q) = means 
that there exists a perfect nesting with the nesting vector 
Q. 

Upon the presence of the particle-hole symmetry, non- 
interacting magnetic susceptibility at Q is given by 
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(4) 



where /(e) is the fermi distribution function. In the 
derivation, we used that D{e)/e is an odd function of 
e so that its principal value integration vanishes. Gen- 
erally, when density of states (DOS) at the fermi level is 
finite, the susceptibility Xo{Q) shows a logarithmic diver- 
gences. On a square lattice, van Hove singularity creates 
a divergence of the DOS, D{e) cx log(e), and hence we 
have a log-square singularity in XoiQ)- 

Using an RPA, the criterion for the AF instability is 
given hy 1 — UxoiQ) =0, where U is the on-site repulsion 
energy. The divergence of XoiQ) implies that an AF 
spin density wave (SDW) occurs at infinitesimal values 
of U. In other words, the critical value for the SDW 
phase transition is Uc = 0. The SDW state is smoothly 
connected to the strong coupling limit U —^ oo, where 
the system is well described by an AF Heisenberg model 
on a bipartite lattice. 

Let us now consider a half-filled Hubbard model on a 
single-layer honeycomb lattice with N.N. hoppings. The 
energy dispersion is given by e±{k) = ±2t'jxy{k), where 
t is the hopping energy. The energy dispersion has a 
particle-hole symmetry due to the bipartite nature of the 



lattice. What is unique in this lattice is that the DOS 
decreases to zero at the fermi level. 



D{e) cx |e| 



(5) 



This may be understood from the fact that, at /i 0, the 
fermi surface shrinks to points at -?C-points (±27r/\/3, tt). 
Using eq. (||) in eq. (4), we see that the singularity in 
XoiQ) is suppressed. Then, RPA gives that the AF in- 
stability occurs at a non-zero value of Ucr-f^ l/xa{Q)- 
The detailed calculation gives Uc/t = 2. 2313) 

In order to examine this result with respect to the 
reliability of the RPA treatment, we perform a quan- 
tum Monte Carlo calculation of the Hubbard model on 
a honeycomb lattice. We anpjy a projection algorithm 
with optimized initial states.liS The ground state is ob- 
tained by I^'g) — exp(— rH)|^'ini) at large enough values 
of r. Here, H is the Hamiltonian of the Hubbard model 
with electron hopping energy t and a repulsive interac- 
tion U. For the initial state |^ini), we use a Hartree-Fock 
solution for the identical lattice with an interaction Umi- 
Note that C/jni (7^ U) can be chosen arbitrarily, and an 
appropriate value for Uini accelerates the convergence to 
the ground state. 

Monte Carlo runs are performed for various system 
sizes with linear dimension L = 6, 8, 10, 12, 14 and 16. 
We typically run 10,000 Monte Carlo steps after 1000 
relaxation steps. AF structure factor is calculated by 



S{Q) 



(6) 



Here N is the number of sites. Staggered magnetization 
Ms is obtained by 



Ms^ = lim S{Q)/N. 



N 



(7) 



Therefore, existence of the AF order can be judged by the 
extrapolation of S{Q)/N to the thermodynamic limit. 

System size dependence of S{Q)/N for various values 
of U is depicted in Fig. 2. We extrapolate the data by 
1/L. Here, finite size effects due to linear spin- wave type 
excitations in the AF ordered states are assumed. From 
the data we see that the extrapolation gives Mg > at 
U/t > 3.7. Meanwhile, at U/t < 3.5 the 1/L extrapola- 
tion gives negative values for Mg^, indicating that 1/L- 
depcndence derived from the linear spin wave theory does 
not hold in this case. From this fact, the critical value £qi- 
the onset of AF order is estimated to be Uc/t ~ S.B.Q 
Thus, quantum fluctuations do not extinguish RPA in- 
stability. They only renormalize the value of Uc/t by a 
factor of about 1.6. As long as the existence of the in- 
stability is concerned, RPA is considered to be reliable. 

Next, we introduce inter-layer hoppings. For a while, 
we assume a half-filled situation for the 2pz bands, of 




Fig. 2. System size dependence of ttie AF structure factor for a 
Hubbard model on a single-layer honeycomb lattice. Lines in the 
figure are least-squares fits for V jt = 3.5, 3.7 and 4.0 from the 
bottom to the top. 



which energy is described by eq. (1). In actual MgB2, 
however, existence of hole fermi surfaces in the 2px^y 
bands leads a deviation from half-filling for the 2pz 
bands. The doping effect will be discussed later. 

As illustrated in Fig. 1, the lattice structure is bi- 
partite. We take intra-layer and inter-layer N.N. hop- 
pings only. An inter-sublattice hopping on a bipartite 
lattice creates a particle-hole symmetry. Indeed, a non- 
interacting electron dispersion in eq. (1) shows 

£+{kx, ky, kz) — £— (fcj;, ky, k^ -\- 7r). (8) 

Note that, since a unit cell contains two boron atoms and 
the first Brillouin zone is already halved by the crystal- 
lographic reason, the AF wavevector in this lattice is 
Q = (0,0, tt). As an inter-layer hopping tz is increased 
from zero, the energy dispersion along the kz axis makes 
the fermi surface expand from the lines in the Brillouin 
zone connecting K and H points, and creates fermi sur- 
faces with finite surface volumes around kz — and 
kz — TT. Therefore, aX tz ^ 0, the DOS at the fermi 
level becomes non-zero. Then, from eq. (4), XoiQ) di- 
verges logarithmically, and the RPA instability occurs 
at infinitesimal values of U. 

This behavior can be understood as a creation of a 
fermi surface nesting by increasing the inter-layer elec- 
tron hoppings. For any fc-points on the fermi surface, 
the relation e+{kx,ky,kz) = —e-{kx,ky,kz -I- tt) = 
holds, which means the perfect interband nesting with 
Q = (0,0, tt). The situation is schematically illustrated 
in Fig. 3. Increase of tz makes the surface volumes of the 
perfectly-nested fermi surfaces increases. 

Enhancement of the nesting behavior by the inter- 



Fig. 3. Schematic illustration of the fermi surfaces and the nest- 
ing vector on a layered honeycomb lattice with finite inter-layer 
hopping terms. The hole (electron) fermi surface around kz = 
{kz = tt) is shown by solid (dashed) gray curves. The nest- 
ing vector Q = (0, 0, tt) is given by an arrow. The dashed line 
represents the Brillouin zone boundary. 



layer hoppings is in a great contrast with usual quasi 
low-dimensional systems where three dimensionality in- 
duces curvatures of the fermi surface and hence decreases 
nesting properties. This unique feature comes from the 
bipartite geometry of a honeycomb lattice. 

Within this band structure, Uc holds strictly only 
at tz = 0, while we suddenly have Uc — Sit tz ^ 0. Uc 
behaves singularly as a function of tz- In other words, for 
a small value of U, SDW order appears as soon as inter- 
layer hopping is allowed. In practice, however, there ex- 
ist particle-hole symmetry breaking terms, e.g. second 
neighbor hoppings as well as a chemical potential shift 
due to carrier doping. Then, there exists a competition 
between the inter-layer hopping term which create the 
fermi surface nesting and the symmetry breaking terms 
which destroys the nesting. As a result, the singular be- 
havior in the tz dependence of Uc will be smeared out 
by the particle-hole symmetry breaking terms, and Uc 
is speculated to remain non-zero. Nevertheless, as long 
as the particle-hole symmetry breaking terms are small, 
we expect that AF instability occurs at relatively small 
values of U. 

If the 2pz bands of MgB2 are half filled, the strong 
nesting behavior of the fermi surfaces creates an AF- 
SDW state even if the Coulomb repulsion is not so large. 
Therefore, it is possible that MgB2 is a nearly- AF metal. 
Self-doped carriers from the 2pj;^y bands to the 2pz bands 
might have destroyed the AF order and make MgB2 
nonmagnetic, being similar with the case for high-Tc 
cuprates in the doped region. At this point, however, it 
is not clear whether such AF is correlated to the mech- 
anism of the superconductivity in MgB2, e.g. whether 
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remnant AF fluctuations assist the phonon mechanism to 
enhance the attractive interaction for the Cooper pair- 
ing. Further detailed studies are required to address this 
question. 

In this sense, lattice structure of MgB2 is special 
because the layer stacking conserves the bipartite na- 
ture. As a counter example, we make a comparison with 
graphite, where stacking of the honeycomb layers are not 
similar to those for MgB2. A carbon atom on an adjacent 
layer lies on top of the center of a hexagon, which is not 
a bipartite structure. Then, the system is stable against 
AF instabilities in the small U region. This explains the 
fact that graphite does not show AF. Coulomb interac- 
tions are expected to be not strong enough to induce a 
magnetically ordered state. 

In order to justify the idea that MgB2 is a nearly- AF 
metal, we propose some crucial tests. From the band cal- 
culation, more precise estimate for the band dispersion, 
either by the tight-binding fits, or by directly calculat- 
ing the fermi surface, will give how strong the nesting 
behavior is in this material. Coulomb repulsion energy 
should also be estimated to clarify whether it is strong 
enough to cause AF states if the 2pz bands are half filled. 
Experimentally, measurement of remnant AF spin fluc- 
tuations in MgB2 is important. A search for AF states 
in various diboride compounds with half-filled 2pz bands 
is also interesting. 

To summarize, ground state magnetism of the Hub- 
bard model on a layered honeycomb lattice at half-filling 
is studied as a model for the B 2pz bands of MgB2. RPA 
as well as MC calculations are performed to study AF- 
SDW singularities. On a single-layer honeycomb lattice, 
AF instability occurs at Uc/t ~ 3.6. On a layered honey- 
comb lattice, inter-layer hopping term induces the fermi- 
surface nesting, contrary to the usual cases where three 
dimensionality destroys nesting behaviors. Within RPA 
scheme, the perfect nesting leads AF to occur at J7c = 0. 
We discuss that MgB2 might be a nearly-AF metal. 

The computations have been performed on the 
massive- parallel personal computer systems of AOY AMA-h 
project ( http : / / www. phys . aoyama. ac . j p / ~ aoyama -I- ) . This 
work is supported by a Mombu-Kagakusho Grant-in- Aid 
for Scientific Research. 
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